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Abstract 

We have measured the relative phase of two Bose-Einstein condensates (BEC) 

using a time-domain separated-oscillatory- field condensate interferometer. A 

single two-photon coupling pulse prepares the double condensate system with 

a well-defined relative phase; at a later time, a second pulse reads out the 

phase difference accumulated between the two condensates. We find that the 

accumulated phase difference reproduces from realization to realization of the 

experiment, even after the individual components have spatially separated 

and their relative center-of-mass motion has damped. 
PACS Numbers: 03.75.Fi, 05.30.Jp, 32.80.Pj, 42.50.Dv 
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The relative quantum phase between two Bose-Einstein condensates is expected to give 
rise to a variety of interesting behaviors, most notably those analogous to the Josephson 
effects in superconductors and superfiuid ^He [|T|. Experiments with condensates realized in 
the dilute alkali gases [0-^ have recently drawn considerable theoretical attention, with a 
number of papers addressing schemes by which to measure the relative phase. Two 
independent condensates are expected to possess |^ (or develop upon measurement [pi|lO|) 
a relative phase which is essentially random in each realization of the experiment. The 
experimental observation at MIT of a spatially uniform interference pattern formed by con- 
densates released from two independent traps llTI is consistent with this view. In this 
Letter, we use an interferometric technique to measure the relative phase (and its sub- 
sequent time-evolution) between two trapped condensates for which the relative phase is 
initially well-defined. This system permits us to characterize the effects of couplings to the 
environment on the coherence |jl2| between the condensates. 

The apparatus and general procedure we use to attain BEC in Rb are identical to those 
of the more recent |TB|,[TB| of our previous two-condensate studies [^^-|T^. We load roughly 
10^ atoms in the \F = l,mi? = —1) (|1)) spin state of ^"^Rb into a time-averaged, orbiting 
potential (TOP) magnetic trap |[TB|. We then magnetically compress and evaporatively 
cool the gas for 30 s, producing a condensate of approximately 5 x 10^ atoms with no 
noticeable non-condensate fraction (> 75% of the gas is in the condensate). The rotating 
magnetic field (i^af = 1800 Hz) is then ramped to 3.4 G and the quadrupole gradient to 
130 G/cm, resulting in a trap with an axial frequency = 59 Hz. The fields are chosen 
to make the hyperfine transition frequency nearly field- independent [0. We create the 
second condensate by applying a short (~ 400 /is) two-photon pulse that transfers 50% of 
the atoms (|-pulse) from the |1) spin state to the \F = 2, nip = 1) (|2)) spin state. The 
coupling drive consists of a microwave photon at 6833.6640 MHz and a radiofrequency (rf) 
photon at 1.0134 MHz; the sum of these frequencies is detuned slightly (~ 100 Hz) from the 
expected transition frequency in our trap PD[. After an evolution time T and an optional 
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second ^-pulse, we release the condensates from the trap, allow them to expand, and image 



either of the two density distributions [15]. The post-expansion images preserve the relative 



positions and gross spatial features of the condensates as they were in the trap p!6| , |2l 



The evolution of the double condensate system, including the coupling drive, is governed 
by a pair of coupled Gross-Pitaevskii equations for condensate amplitudes $i and $2: 

= {T + V, + U, + Uu) $1 + (1) 

and 

th^ = {T + V2 + Vms + U2 + U21) $2 + ^e*--*$i (2) 

where T = — (/;,^/2m)V^ is the kinetic energy, m is the mass of the Rb atom, Vms is the mag- 
netic field- dependent hyperfine splitting between the two states in the absence of interactions, 
condensate mean-field potentials are Ui = A-nh'^aini/m and Uij = A-nh^aijUj/m^ rii = 
is the condensate density, and the intraspecies and interspecies scattering lengths |l5|Jl6 



are and aij = aji. For the trap parameters given above, the harmonic magnetic trap- 
ping potentials Vi and V2 are displaced from one another by 0.4 fim along the axis of the 
trap [|l^. The coupling drive is represented here in the rotating wave approximation and is 
characterized by the sum of the microwave and rf frequencies, ujj.f, and by an effective Rabi 
frequency where 

{27r ■ 625 Hz, coupling drive on; 
(3) 
0, coupling drive off. 

Phase-sensitive population transfer between the |1) and |2) states occurs with the drive on, 

but the two condensates become completely distinguishable once the drive is switched off. 

The first |-pulse [Fig. |l](b)] creates the [2) condensate with a repeatable and well-defined 

relative phase with respect to the [1) condensate at t = 0. The relative phase between 

the two condensates subsequently evolves at a rate proportional to the local difference in 

chemical potentials between the two condensates 0^21 (r,t), which in general is a function 

of both time and space. Couplings to the environment |2^ can induce an additional (and 
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FIG. 1. A schematic of the condensate interferometer, (a) The experiment begins with all of 
the atoms in condensate |1) at steady-state, (b) After the first ^-pulse, the condensate has been 
split into two components with a well-defined initial relative phase, (c) The components begin to 
separate in a complicated fashion due to mutual repulsion as well as a 0.4 //m vertical offset in the 
confining potentials (see also Fig. 3 of Ref. |jl^). (d) The relative motion between the components 
eventually damps with the clouds mutually offset but with some residual overlap. Relative phase 
continues to accumulate between the condensates until (e) at time T a second ^-pulse remixes the 
components; the two possible paths by which the condensate can arrive in one of the two states 
in the hatched regions interfere, (f) The cloud is released immediately after the second pulse and 
allowed to expand for imaging. In the case shown, the relative phase between the two states at 
the time of the second pulse was such as to lead to destructive interference in the |1) state and a 
corresponding constructive interference in the |2) state. 
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uncharacterized) precession of the relative phase, leading to an rms uncertainty in its value 
A(y9diff [^,0. After an evolution time T, therefore, the condensates have accumulated a 
relative phase Jq uj2i{r, t) dt + A(y9difr(T). During the same time, the coupling drive accumu- 
lates a phase UrfT. A second |-pulse [Fig. 0(e)] then recombines the |1) and |2) condensates, 
comparing the relative phase accumulated by the condensates to the phase accumulated by 
the coupling drive. The resulting phase-dependent beat note is manifested in a difference 
in the condensate density between the two states. Immediately after the second pulse the 
density in the |2) state (n2/) is 



'^2/(r) = ^ni{r) + ^n2{r) + \jni{f)n2{r) cos U uj2i{r,t) dt \ - u^^T + Alp^^s{T) 



■ (4) 



In this equation, rii denote the densities prior to the application of the second |-pulse. 
The interference term in Eq. ^ shows that measurement of n2/(r) in the overlap region is 
sensitive to the relative phase. Each realization of the experiment (with a freshly-prepared 
condensate) yields a measurement of the relative phase for a particular T; by varying T, we 
can measure the evolution of the relative phase. 

At short times T, for which the overlap between the condensates remains high, varying 
the moment at which the second |-pulse is applied causes an oscillation of the total resulting 
number of atoms in the |2) state. The oscillation occurs at the detuning frequency 5 = UJ21 — 
ujj-i and is completely analogous to that observed in separated-oscillatory-field measurements 
in thermal atomic beams p8[ or in cold (but noncondensed) atoms in a magnetic trap p9 



The fringe contrast, initially 100%, decreases as the condensates separate. After ~ 45 ms 
the relative center-of-mass motion damps and comes to equilibrium, leaving the components 
with a well-defined overlap region at their boundary, as shown in Figs. |l](d) and 0(a); see 
also Fig. 5(b) of Ref. [|I6|. Application of a second |-pulse at T > 45 ms results in a density 
profile in which the interference occurs only in the overlap region; see Figs. |l|(f) and|^(b). 
We look at the density of atoms in the |2) state at the center of the overlap region |^ in 



order to examine the intriguing issue of the reproducibility of the relative phase accumulated 
by the condensates during the complicated approach to equilibrium. If the phase diffusion 
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FIG. 2. (a) The post-expansion density profiles of the condensates in the steady-state attained 
after a single |-pulse. These density profiles vary little from shot-to-shot (and day-to-day), (b) 
The density profiles after the second |-pulse. The density in the overlap region depends on the 
relative phase between the two condensates at the time of the pulse; in the case shown, we observe 
constructive interference in the |2) state and destructive interference in The patterns in (b) 
are much less stable than those in (a), possibly as a result of unresolved higher-order condensate 
excitations, issues associated with the expansion, or technical instabilities of the apparatus. 
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term in Eq. ^ is so large that the uncertainty is greater than vr, then repeated measurements 
for the same values of T will yield an incoherent [i.e., random) ensemble of interference 
patterns. In the opposite extreme, {i.e., very little phase diffusion), repeated measurements 
will give essentially the same interference pattern at T in each experimental run. We plot the 
optical density in the center of the overlap region as a function of T in Fig. |^, and observe an 
oscillation at the detuning frequency with a visibility of approximately 50%, corresponding 
to an rms phase diffusion A(pdis{T) ^ |. At longer times the maximum contrast observed 
in a single realization of the experiment decreases slightly, possibly due to the increasing 
presence of thermal atoms as the condensates decay. 

The stable interference patterns show that the condensates retain a clear memory of their 
initial relative phase despite the complicated rearrangement dynamics of the two states fol- 
lowing the first |-pulse. This is rather surprising, since the center-of-mass motion of the 
double condensate system is strongly (and completely) damped, and, in general, decoher- 
ence times in entangled states tend to be much shorter than damping times pT|-p3|. The 
intuition one develops in understanding few-particle quantum mechanics may not apply to 
experiments involving condensates. The phase between the two condensates seems to pos- 
sess a robustness which preserves coherence in the face of the "phase-diffusing" couplings to 
the environment. 

We have read out the relative phase of two Bose-Einstein condensates using a time- 
domain version of the method of separated oscillatory fields. We observe the persistence of 
phase memory in this "condensate interferometer," despite the presence of damping and the 
complicated rearrangement of the two condensate components. We have established that 
the time scale for phase diffusion in this system can be longer than 100 ms. The double- 
condensate methods we have developed will be applicable to other experiments which explore 
phase diffusion as a function of condensate parameters including temperature, number of 



atoms p^-p7|, and collision rates |34|. Collapses and revivals of the "memory" of the 



relative phase are predicted (jTU],^ at time scales which may be experimentally accessible 
should environmentally-induced diffusion effects remain small. Our methods will also allow 
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FIG. 3. The value of the condensate density in the |2) state is extracted at the center of the 
overlap region (inset) and plotted (a) as a function of T. Each point represents the average of 
6 separate realizations and the thin bars denote the rms scatter in the measured interference for 
an individual realization. The thick lines are sinusoidal fits to the data, from which we extract the 
angular frequency UJ21 — uj^f. In (b), the frequency of the coupling drive Wrf has been increased by 
27r ■ 150 Hz, leading to the expected reduction in fringe spacing. 
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us to examine other phase-related phenomena, such as phase-locking and analogues of the 
superconducting Josephson junctions [p5| . 

We gratefully acknowledge useful conversations with A. J. Leggett, as well as with the 
other members of the JILA BEC collaboration. This work is supported by the ONR, NSF, 
and NIST. 
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